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ABSTRACT 


This report examines the effect of a series of evenly spaced moving 
loads on an elastic rectangular plate that is clamped at all edges. Two types 
of loads are treated: (1) a uniformly distributed pressure over a fractional 
length of the plate and ( 2) an impulsive load. All loads across a plate are con- 
sidered uniforpa. 

The solution of the partial differential equation of vibration of a plate is 
assumed in the form of a double series with the generalized coordinates solved 
by using the Laplace transform method. Viscous damping is included. 

Based on a two-term approximation, the steady-state dynamic response 
of the plate is obtained in analytical form from which the upper bounds of the 
maximum deflection and maximum bending stress are formulated. Numerical 
examples are given to illustrate the effect of the thickness of the plate, the 
aspect ratio, and the velocity of loads on the dynamic response. Three types 
of resonance conditions are derived. 
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DEFINITION OF SYMBOLS 
Definition 
Length of plate 

Coefficient of expansion, equation ( 8) 
Generalized coordinates, equation (6) 
Width of plate 

Coefficient of expansion, equation ( 8) 
Velocity of loads 

Nondimensional velocity of sound 

Equation (5) 

Eh^ 

12 (1-v 2 ) — * lending stiffness of plate 
Modulus of elasticity 
Thickness of plate 

Coefficient of expansion, equation (11) 
Length of distributed loads 
Laplace transform operators 
Intensity of load 
Dynamic load function 
a / b, plate aspect ratio 
Laplace transform variable of t 
Time variable 
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DEFINITION OF SYMBOLS (Concluded) 


Symbol 

u(t) 

w(x, y, t) 

x, y 


a . 

1 


n 

6(t) 

£ 

0 

n 

A 

n 

M 

n 


(z) 
^ n 


co 


n 


9825, oi 2 = 1.0008) 


Definition 
Unit step function 
Dynamic deflection of plate 
Coordinates of plate 
Constant of beam function ( a i = 

=A 2 / Aj , amplitude ratio 

- \ A * 

Dirac function 

Damping factor 

= c A 

n 

Eigenvalues of a clamped-clamped beam (Ai=4. 7300=1. 50567T 

A 2 =7. 8532 =2.4997t r ) 

= 0 / co 

n n 

Poisson’s ratio 
Mass density of plate 
Bending stress 

Time interval between two consecutive loads 
Eigenfunction of a clamped-clamped beam 
= co /T, natural frequencies of plate 




D 


Note: Barred symbols denote dimensional quantity, and unbarred symbols de- 
note nondimensional quantity. 


vi 
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ON DYNAMIC RESPONSE OF A RECTANGULAR PLATE 
TO A SERIES OF MOVING LOADS 

SUMMARY 

This report examines the effect of a series of evenly-spaced moving 
loads on an elastic rectangular plate that is clamped at all edges. Two types 
of loads are treated: (1) a uniformly distributed pressure over a fractional 
length of the plate and (2) an impulsive load. All loads across a plate are con- 
sidered uniform. 

The solution of the partial differential equation of vibration of a plate is 
assumed in the form of a double series with the generalized coordinates solved 
by using the Laplace transform method. Viscous damping is included. 

Based on a two-term approximation, the steady-state dynamic response 
of the plate is obtained in analytical form from which the upper bounds of the 
maximum deflection and maximum bending stress are formulated. Numerical 
examples are given to illustrate the effect of the thickness of the plate, the as- 
pect ratio, and the velocity of loads on the dynamic response. Three types of 
resonance conditions are derived. 


INTRODUCTION 


The primary objective of this report is to determine the dynamic re- 
sponse of an elastic rectangular plate subjected to a series of evenly-spaced 
moving loads. Both a uniformly distributed pressure over a fractional length 
of plate and an impulsive load are treated herein. In fact, an impulsive load 
is a limiting case of a uniformly distributed pressure, when the length of the 
distributed pressure becomes infinitesimal while the product of length and 
intensity takes the finite value, P 0 . 

There are some practical applications of the preceding mathematical 
model to acoustic problems or stress analysis of a shell structure in an aero- 
dynamic flow. Consider a plate that is an idealized panel of an airplane wing, 
or a skin panel of a large shell structure, while the moving loads resemble 



shock waves or pressure disturbance originating from a noise source. If the 
loads travel at very high speeds relative to the plate at sonic velocity for in- 
stance it takes less than one hundredth of a second for the load to go across a 
plate. When this is in the same order of magnitude as the period of vibration, 
resonance may take place under certain relationships between velocity of the 
loads and the natural frequency of the plate. Furthermore, conditions of res- 
onance may be related to both velocity and frequency of the moving loads. A 
determination of the dynamic response of the plate at resonance conditions can 
also be made if the viscous damping factor is known. To a structural designer, 
all this information may be helpful for the determination of the optimum aspect 
ratio and thickness of the panel. 

The assumption is made that the loads travel in one direction only and 
that the space between two consecutive loads is greater than the length of plate. 
Otherwise, the principle of superposition may apply. The coordinates 
and dimensions of the plate and loads are shown in Figure 1 . Analysis begins 
with the vibration of a rectangular plate with clamped edges. 



FIGURE 1. COORDINATES OF PLATE 
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EQUATION OF MOTION 


The response of an elastic plate, w(x, y, t), to a dynamic load, p(x, y, t) , 
is governed by the equation of vibration of a plate 


D(w + 2w 1- w ) +phw-- = 

xxxx xxyy yyyy tt 


P(x,y,t) 


( 1 ) 


— 3 w — — 

with boundary conditions w = — = 0 for clamped edges. The subscripts x, y, 

and t denote the partial derivatives of w. A symbol with a bar represents a 

dimensional quantity; the unbarred symbols denote the nondimensional quantity. 

The function p(x, y, t) of the two types of moving loads to be treated can 
be expressed 

For impulse loads (L— ■ 0) 

OO 

p(x,y,t) = P 0 Yj 6 ( t-x/c-kr) , (2a) 

k=0 


For distributed loads with length L ( L < a) 

OO — - 

p(x,y,t)=P„ Y [u( t-x/c-kr)- u(t- ^4 -kr)] (2b) 

k=0 C 

where <5 (t) is the Dirac function and u(t) is the unit step function. 

For convenience of manipulation, rewrite equation (1) in the nondimen- 
sional form 


4- (w +2r 2 w + r*w ) +w,, = p(x,y, t) 

X? xxxx xxyy yyyy tt 


(3) 


The barred and unbarred symbols are related as follows 
w=w/a, x = x/a, y= ry/a, r= a/b L= L/a 


A / A 2 \ 3 

T= n/ 4)h/mf , t= t/T, c=cT/a, r =t/T, P 0 = P 0 T 2 Aph = P& - 

The definition of is investigated in the succeeding section. 
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METHOD OF SOLUTION 


Observing that p(x, y, t) has a simple expression in Laplace transform, as 
a first step toward the solution of equation (3) , we apply Laplace transformation 
to the partial differential equation with respect to the time variable. It follows 
immediately that 


rV(W + 

X* xxxx 


2r 2 W + 
xxyy 


r 4 W 


yyyy 


) + s 2 W = P 0 C (s)e 


-sx/ c 


(4) 


where 


W = W( x, y, s) = L{w(x, y, t)} 

-TS) 


C(s) = 


[1/ (i -e 
I < 1 — e -slj/ 0 


for impulse loads 


-TS 


)/ s(l - e ) for distributed loads 


(5) 


Note that the initial conditions w (x, y, 0 + ) = w (x, y, 0 + ) = 0 have been taken. 

t 

A well-known approach to the solution of equation (4) is to assume the 
solution in the form [ 2 ] 

oo oo 

W (x,y, s) = ]] Tj A nm (s)< ^n (x) V (y) (6) 

n=i m=l 


The term 0 .( z) is the eigenfunction of the differential equation of a clamped- 
clamped beam. 


d 4 0 

d? 


-X 4 0 = 0, 


with the boundary conditions 0(0) = 0 (1) =0’ (0) =0' (1) =0. This assumed 
solution satisfies the boundary conditions along all four edges but does not 
satisfy the differential equation. However, the differential equation can be 
satisfied approximately by using the generalized Fourier series expansion. 

The sequence 

0 . ( z) = coshX .z - cosX.z — a . (sinhX.z - sinX.z) i = 1,2,... (7) 

r i l ill i 

a . = (coshX .- cos X.)/ (sinhX. - sinX.) 

l i l i l 


4 


f 


forms an orthogonal, complete set in the interval 0 to 1 . This permits expansion 
of the term on the right-hand side of equation (4) in terms of 4> ^(x) <p m (y) 

e - sx/o = £ a n (s)b m ^(x)^ m (y) (8) 

m, n=l 


where, as shown in Appendix A, 


a (s) 
n 


sUe 2 

n 


20 ^ 


n , J% n, _ v -s / c 

— -(-1) (s + 0 )e 

s + 0 n 

n 


, 0 = cX 
J n n 


b 


4a / X m odd 
m m 


^0 m even 

To take into account the damping effect on the vibration, assume that the damping 
is viscous and has a viscous damping coefficient of 


e 

nm 


<P n (x)<p 


m 


(y). 


Substituting from equations ( 6) and ( 8) into equation ( 4) and adding the 
damping term leads to 


2 

n, m= 1 


A (6 4 +r 4 6 4 +e s + s 2 )<p (x)<p (y) + 2r 2 6 2 fi 2 <p " (x) <p " (y) 
nm n m nm ^n ma n mm ' ma J 


= P 0 C(s) £ 


n, m=l 


a (s)b <p ( x)<p (y) 
n m Y n ^ m J 


(9) 


where 


<?<t> (z) 

<*>[' (Z) = 1? d? and 6 i = V Xl 

i 


It is now necessary to expand <p ^ ( z) in terms of <p . ( z) ; this gives 

oo 

<j»;'(z)= 2 k * (*) 

J i=l J 


( 10 ) 
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The coefficients of expansion are 


[4] 


k..= j <t>\' (z)0.(z)dz = 4[l+(-l) 1+ ^ X 2 . (a.X. a.\.)/(\f -\?) i * j 

ij J o J 1 L J 1 1 1 " J J 1 J 


k.. = f (f>."( z) 0 .(z)dz = (2 - o;.X.)o!. A. 

11 J n 1 1 llll 


( 11 ) 


Substitution from equation (10) into equation (9) results in a set of an infi- 
nite number of simultaneous algebraic equations in the generalized coordinates 


A (s) 
nm 


OO 

(< 5 4 + r 4 5 4 +e s+s 2 )A + V 21^6? 6? k .k .A.. = P<>a (s)b C (s) 
n m nm nm . . ^ ^ i j m mj ij w n m 


( 12 ) 


n, m = 1,2, 


Evidently, if a large number of equations is used to obtain the solution of 
equation (12), the expression of A nrn (s) becomes combersome, and the task of 
finding inverse Laplace transforms becomes prohibitive. To serve the purpose 
of the present study, a simple approximation of m = 1 and n = 1 and 2 is made. 
The solution of A^ (s) and A 21 (s)can then be readily obtained from the reduced 
system of equations 


1+ r 4 + 2r 2 k 2 11 + e^ s + s 2 0 

0 6f + r 4 + 2r 2 ^i k 22 kn +e 22 s+s 2 



A il 


a t (s) 



=P 0 biC(s) 



A 2 i 


a 2 (s) 


Notice that equation (11) gives k 12 = k 21 = 0; consequently, the first two natu- 
ral frequencies of the plate can be written out directly from equation ( 13) . 


(13) 


oj — <j0 /T — 2. 267o> 


n 


n 


(—si 

1 D ) 

la 2 ^ 

Ph / 


00 


= \ / 6 4 + r 4 + 2 6 2 i^kiik 
n n n 11 i 


(14) 


nn 


The nondimensional frequencies versus the plate aspect ratio given by equation 
(14) are plotted in Figure 2. 
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2 3 4 

r = a/b 


. NONDIMENSIONAL FREQUENCIES 
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To further simplify the solution of A lt and A 21 , assume that the damping 
factors for the first and second modes are equal. 


% = i e ii wi = 2 e 22^2 

Hence, the damped natural frequencies of the plate become 


CO =co n/ 1 -C, 2 
nd n 


Express the solution of equation (13) in the form 

-s/ c' 


where 


A nl (s) =P 0 b 1 (c/co 2 n ) [R n (s) +Q n (s)e' s/C jc(s), 

29 2 

R n <s) = 


Q n (s)--(-l) 


n (s +Q n ) 

[ (S+ ^ n )2 ^ 2 nd] (s2 + 0 n ) 


( 15 ) 


( 16 ) 


(17) 

( 18 ) 


Determination of the inverse Laplace transform for the two types of loads will 
be given separately. 


Solution for Impulse Loads 


By using the method of partial fractions, R n (s) equation (17) and Q n ( s ) 
of equation (18) can be written as a sum of five and four functions, respectively 


V s * = i R „k F „k (s > • V s) =|v F nk< s) ■ 


(19) 


The constants, R , , Q , , and the function , F . (s) , are shown in 
nk nk nk 

Appendix B. For impulse loads, C(s) of equation (16) is as noted in equation 
(5), 


C(s) = 


1_ 

1 - e 


-ts 
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It has been shown [5] that the inverse Laplace transform of a periodic function 
may be expressed in the form: 


-l 

L 


y(s) |= M'S) - * (t) , 
l-r™J 


Kt < t < (K+i) r, K = i, 2, ... 

A 

- T < t < 0 


( 20 ) 


in which ipd) is a function of a discrete variable^, and ^(t) is the same func- 
tion of the continuous variable t. The Laplace transform pairs of F and ip of 
equation (20) are given in Appendix B. As derived 


± L- 
c 


i ,-l 


V s) 


1 - e 


-rs 


= R 


no 


i p (fy - ^ (t) 

no no 


+ f (t) - f (t) + g ft) - g (t) 
nr nr nr nr 


Q„(s) 


1 - e 


- TS 


■= f (fy - f ( t) + g (fy - g ( t) 

nq nq b nq to nq 


(21a) 

(21b) 


Finally, the inverse Laplace transform of equation (16) may be writ- 
ten in the form 


"I- 


L_1 Kil (s) } =P ° bl (c/c °n ) K 


(t) - f n (t) +g n (t) - g n (t)] , 


(22) 


where 


f (t) = H $ (t) +f (t) -(-l) n f (t--^-) u(t- -7) 

n no no nr nq c c 


(22a) 


g (t) = g (t) - (-i) n g (t- — ) u (t- — ) 
n nr nq c c 

A A A A n A 1 

f n (t) =R „nt (t) + *<*> " f (t + LT — ~ ) 

n no no nr nq c 

A /A\ / a \ ^ /A ^ \ 

g n (t)=g „r (t) -<- 1) *nq (t + lT ‘T ) 


^ 1 1 
when - t < t < - r+ — 

c 


when - r + < t < 0 

c 


L = 1 


(22b) 

(22c) 

(22d) 


t = 0; and g (ty = g $) , £ $) =f $) 
n n n n 
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Substitution from equations (B-ic) and (B-id) into equation (22b) yields 
g (t) = \ J 2B [1- (-l) n sinA ] csc^Q t sin[© (t+|r)+y+T? ] for t>— (22b') 


n 


n 


'n n 


Fora steady-state solution, f n (t) and g^(t) can be disregarded because 

the former decreases exponentially with time while the coefficient of the latter 
is of a highe 

5© t = k7T . 


is of a higher order of smallness (sinA = (-l) n ) except in the neighborhood of 

n 


n 


As shown in Appendix B equations (22c) and (22d) can be rewritten as follows 

(23a) 


«. A 

AA A A- fct)t A 

f(t)=R 0 ^ 0 (t) +Fe COS ( CO t- <p- ip) 


M =J 


^2Bsin [©(^+ ir) +y + T))] when -t<^<-t+-~ (23b) 

0 when - r + — < f < 0, except in the neighborhood of 
c 

J0.T = k7T 


n 

Hence, the two-term steady-state solution of the dynamic response of the plate is 

2 

S 

n= 1 


z 

w(x,y,t) = Pobj (c/co n 2 ) [^ n (^ + S n ^)] <£ n (x)0i(y) 


(24) 


The Upper Bound of the Maximum Deflection and the Maximum Pending Stress. - 
The maximum deflection and the maximum bending stress of the plate can be ob- 
tained by using equation (24) , provided the greatest absolute value of f({Vt-g $) 
is determined. It is rather tedious to determine the greatest value of £ (fy (^) ; 
however, the upper bound of this value, as denoted by f * + g* , can be written 
out immediately 


f* +g* = J 


Roi/'o (" t) + fte ScoT + n/2B 


A 1 

- T < t < - T + 

C 


Ro^o 


, 1 . A , £cv(t- ) 1 A . 

(-T + — )+Fe c -r + <t < 0 


(25) 


Hence, the upper bound of the maximum deflection may be expressed in the form 


w* s aPobiAj^ i(x* ) + p<p 2 (x^* )] <p t (|) = WoAt 


(26-1) 


10 



where 



* = <M** ) +P <M** ) (28) 

w w 

The location of the maximum deflection which is denoted bvx* in the above 

w 

equations can be determined numerically as a function of 0 from the equation 

I <p\(x) I + 62/3 I02 1 (x) |=0 (29) 

1 n 

where <p ' (x) = - — — — . The graphs of x* vs 0 and # vs 0 are shown in 

n X dx w 

n 

Figure 3. 



0 0.2 0.4 0.6 0.8 1 0 1.2 

0 


FIGURE 3. $ and x * VERSUS 0 

w 
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The upper bound of the maximum bending stress in x-direction can be 
determined by using the stress equation 

_ Eh [ 8? w O 2 w 

°x - 2(l-i/) + v dy 2 _ 

with w(x, y, t) given by equation (24) . Analogous to equation (26-i) , express 
the upper bound of bending stress ct * in the form 


5* g <7o Aj [s±(x^ ) + p Szix* )J 4>i (i)= (To AiS (30-1) 

where 

(7 0 =P 0 biX 2 i 0i(i) Eh/ 2 (i-v )a = 0. 47i8P 0 a 2 / h 2 (forv=|-) 

s n ( x )= |fi*0£( x ) + vr2 ^ n ( x > 0 1' C i ) / 0 1 (i ) t ( 31 ) 


S = Si(x*) + p Sa (x*) (32) 

(T <7 

The location of the maximum bending stress x* can be determined numerically 
as a function of p from the equation 


dSj(x) 

dx 


+ 


/3 


dS2(x) 

dx 


= 0 


(33) 


or 

4> i"(x) - 0.76545 i/r 2 <^'(x) + 4. 5768/3 [^ 2 '" (x) - 0. 2777 i/i 2 0 ’ 2 (x)] = 0 

i rf *n 

where <f>'"(x)= — -g '- j- . The graphs of x * vs /3 and S vs p are plotted 

n a cbr a 

n 

in Figure 4. 

To summarize, the procedure of computing the upper bound of the max- 
imum deflection and the maximum bending stress is: 

i. Determine from equation (14) or Figure 2 
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FIGURE 4. S AND x* VERSUS p 

cr 

2. Compute f* and g* by using formulas in Appendix B and equation (25) 

3. Obtain A and B from equation (27) 

n 

4. Read $ and S from Figures 3 and 4, respectively 

5. Calculate w* and a* from equations (26-1) and (30-1), respectively 

X 

Three Special Cases, u » 1, u « 1 and / u =1.- The three extreme cases, 

n n n 

namely, c»co/\ , c«c o /\ andc=u> /\ , deserve special attention, 
n n n n n n 

Physically speaking, the first two conditions say that the frequencies of the ve- 
locity waves produced by the moving loads are very high and very low, respec- 
tively, in comparison with the natural frequency of the plate, and the third is a 


A •} 


AO 



condition of resonance at which the magnitude of the dynamic response has the 
greatest value. For these three cases, the constants given by the table in Ap- 
pendix B and equations (21) and (B-2) have the approximate expressions given 
in Table I. 


TABLE I. APPROXIMATIONS FOR THREE SPECIAL CASES 



fX » 1 

n 

IX « i 

n 

M = 1 
n 

B 

i/M 4 

1 

1/4 f 2 

C 

i/m 4 

1 

1/ 8t, 2 

D 

-M 2 

1 

2£ 

E 

V 

1 

-2C 

F 

M 3 

3f, -M 

2^ 

G 

A* 3 

?. - M 

24 

Ro 

1/m 2 

1 

i/2(l-U 

R 

i/()M 

m 2 /() 

V2/ 4£,(l-t,) ( ) 

Q 

l/2()/* 

1/2 () 

V2/8£, ( ) 


2 71- 

3^ -ix 

4^r 

e 


TT+^-IX 

4?r 

A 

n 

1 4 


c 

A 

[2-( -1) ]Q 

Qe -t«(ar-3) 

Q* ; 

A 

D 

2sin 2 ^u> d r 

1 

D 


Note that in the above table ( ) = cosh t, lot - costo and 


Q*ssQ - 




2 -2 £co( n -~TT) , i „ n, 4 v - £w(it - ~) 

+ e c -(-1) ( ~ ~g ) e c cosco d (tr- 


Recall the original assumption that £, is much smaller than one. 
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O | M* 



The values of co /X for plate aspect ratio r = 1 to 4 are as follows 
n n 


r = 1 r = 2 r = 3 r = 4 

n = 1 0.341 0.933 1.975 3.450 

n = 2 0.419 0.727 1.326 2.205 

Case I fi »1 orc»w /X 
n n n 

Since V 2B + R 0 is negligible in comparison with f*e ^ WT in equation (25) , 
reduce equation (27) to 


A n“L 2 - ( ' 1>n ] 


.2 1 


X co 1 ( cosh £ co r - cosco , t) 
n n n nd 


(34) 


This expression shows that the magnitude of the dynamic response is independent 
of the velocity of the moving loads when c » co / 

Case II fx < < 1 or c < < co /X 
n n n 

Since t >1 / c or r is large, use the approximation 


Scot 


/ 2 ( coshSco t - cos co ,t) = 1 

d 


This results in 


A s 
n 


1.W2 + e 


-Sco (r- s')" 1 


c/co 


n 


(35) 


Thus, if the third term in the above equation is much smaller than one, the 
magnitude of the dynamic response is directly proportional to the velocity of 
the moving loads when c < < co^/ 


Case III a* = 1 or c = co / X . 

n n n 

The magnitude of the dynamic response at resonance can be written in the 


form 

c/co 2 (36) 

n 

At resonance, A is approximately inversely proportional to the damping factor £. 


A = 
n 


2 (1 - £) 


n /~2 

2£ 


D Q* 
n n 


SCO T 
e n 


1 



Solution for Distributed Loads 


The inverse Laplace transform of equation (16) for the case of distributed 
loads becomes more cumbersome than for the impulse loads because of the com- 
plexity of C(s) given by equation (5) . It will simplify the problem somewhat if 
the damping factor is assumed equal to zero. Now, rewrite equation (16) in the 
form 


A^s) = 


Pnbj 

X CO 2 

n n 


R (s) -(-l) n Q (s)e' 
n n 


s 

c 


-sL/c 


1 - e 


-TS 


where 


AR(s) = 


20 3 c o 2 


s( s 2 + w 2 ) ( s 2 + e 2 ) ( s + e ) 


(37) 


_2_ 2j a 3 ( jxs + co \ 1 /s + 9 \ 

s + co 4 (l -M 4 ) \s 2 + co 2 / 1 - M 2 \s 2 +9 2 / 


XQ(s) = 


0o) 2 (s + e) 
s( s 2 + CO 2 ) ( s 2 + 0 2 ) 


1 

1 +M 2 



(38a) 


co 2 (l -M 2 ) 


MS - co 
s 2 + co 2 



(38b) 


It follows immediately that 

f 

r(t) = L _1 {R( s)} = 2 +Rsinco(t + y) +R' sin© (t + 4»r ) +R"e (39a) 

q(t) = L _1 {Q(s)} = 1 + Qsinco (t - ■>) + Q'sin0(t - ) (39b)) 

where 


R = 


2 ix \} 1 + m 2 

4 4 

CO ( 1 - M ) 


Q = 


M *s/ 1 + M 2 
co 2 (1 - m 2 ) 


- R' = Q' = \T2/(1 - m 2 ) R” = -l/(l + M 2 ) y = tan -1 M 
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Knowing the inverse Laplace transform of R ( s) and Q ( s) , the inverse 

n n 

Laplace transform function of equation (37) denoted by 
A n (t) = L-KA nl ( S ) } 

can be illustrated diagrammatically (Fig. 5) . This function can be expressed 




FIGURE 5. TIME HISTORY OF r(t) AND q(t) 

as sectional, continuous functions in four regions, a to d . With the aid of 
Figure 5 and the shorthand notation 

K 

f(K, a) = Yj f(t-kT-a) u(t-kr-o;) 
k= 0 
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write A ( t) as follows 
n 


Region a. Kr<t<Kr+ L/ c 
Pnb 


A n (t) = > 

n a 


BV lr 
. CO n 

n n 


( K, 0) - r ( K - i,“ ) - ( -1) n [q (K-l,~) 
n c n c 




( 40 ) 


Region b. Kr + — < t < Kr + — 
c c 


A (t) = ^ nb 
n X 


3 0 b 1 

“2~ r 

or n 

n n 


(K, 0) -r (K,-~- ) -(-l) n [q (K-i, -7) 
n c n c 


-q (K-i, ] 

n c 

Region c. Kr + — < t < Kr+ * 1 — 

c c 


A (t) = 
n A 


- (r (K, 

N. or n 

n n v 


( K, 0) - r (K, -7-) - ( -i) n [q ( K, - — ) 

II C II c 


(40b) 


( K-i, )]} (40c) 

n c I 

Region d. Kr+ 1 + - <t< (K + i) t 
c 

A (t) = (r ( K, 0) -r (K, 7- ) - (-i) n [q (K, -7) 

n Aeo n nc nc 

n n 1 

-q^K, ^-^-)]J (40d) 

Making use of the formulas derived in Appendix C, combine the preceding four 
equations into a single equation of variable t 0 which is defined 

t 0 = t - Kr with K = i, 2 , . . . and K t< t < ( K + 1) r 
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A(t)= -^r 

Xlo 


2snr^- [esc |ojr^Rsin|Mo; r cos{ o> [ t 0 + (K- ~^) t- -^-]+ t) 


- (-i) n QsiniPcoT cos{o; [t 0 + (K- r - ] -y}) 

+ Rcos[u> (t 0 + Kr- +y] 

- (-l) n Qcos[o;(t 0 + Kr-^^-) - y]] 

+ 2sin^XL[csc2cAT R' ^singMeXr cos{ cX [t 0 +(K- * ) r - ~ ]+ jir } 

- ( -1) n Q’siniPcXr cos {cX [t 0 + (K- | tt }) 

+ R’cos [cX (t 0 + Kt - ~~ ) + iir] 

- ( -1) n Q'cos [cX(t 0 + Kr - - a tt ]] 

+ R"e - c M t o + Cr) (1 _ e XL ) fl _ e “ c MM+ 1)^ ^ _ e "CXr 

+ C[ 2 +Rsin(cot 0 +y ) + R'sin (cXt 0 + \ n) + R"e c ^ to ] 

- ( -1) n D [l +Qsin [co (to - ~) - y ] + Q'sin [c X(to - ~ ) - 3 ^ l] J (41) 
The values of the constants M, P, C, and D in the above equation are as follows: 


k 

M 

P 

C 

D 

a 

O 

> 

VI 

o 

4-» 

VI 

o 

K-l 

K-l 

1 

0 

b 

O 

VI 

4? 

VI 

o 

K 

K-l 

0 

0 

c 

l/c < t 0 < ( L + l)/c 

K 

K-l 

0 

1 

d 

( L + l)/c < t 0 ST 

K 

K 

0 

0 


Note that for brevity all the subscripts "n" are omitted from the symbols A(t) , 
R, R', R" Q, Q’, X, yandco. 

For the case K < 1, this problem can be treated as a plate subjected to a 
single moving load and is given in the next section. 


19 













( 42 ) 


Now, the dynamic response of the plate is 


2 

w(x, y, t) = aPflbj £ A^(t) <p ^(x) <p t (y) 
n=l 


And 


B n (x,y ) [6 n ^n (x )^ l( ^ + v ^ <P n ( x ) <P 1 (y)) ( 43 > 

Then, the dynamic bending stress in x-direction is 


cr (x, y, t) 

X 


PnbjEh 
2(1- v ) a 


2 


Z 

n=l 


n 


(t) B n (x,y) 


(44) 


3jC 

Using the symbol A to represent the maximum absolute value of A (t) and 

# s|c ^ ^ 

/3 = A 2 /A t , obtain the upper bound of the maximum deflection 

w* = WqAj 4> (26-2) 


and the upper bound of the maximum bending stress 
(7*= CTqA* S 


(30-2) 


where the symbols w 0 , <j 0 , 4> and S are defined previously in equations (26-1) 
and (30-1). 

It is difficult to tell at what time region A* occurs. This suggests a pro- 

n ’■ 

cedure of computing the greatest value of A^( t) in each region and choosing the 

largest one among them. Let F to F denote the greatest value of A (t) in 

a q n 

regions a to d; hence 

A= :< = Max (F , F , F , F ) (45) 

n abed 


A simple but crude method of determining F is to let all the sines and cosines 
of equation (41) equal to one and add the absolute value of all terms. This can 
be refined by first combining the cosines of the same frequency and then summing 
the absolute values term by term. 
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A Single Moving Distributed Load - For a single moving distributed load, the 
load function given by equation ( 5 ) becomes 

L {p (x, y, t) j = — e (1 - e ) ( 46 ) 

s 

Consequently 

L _I { A (s) } = . Pflbl 2 1 I r (t) - r (t-L/c)u(t-L/ c) 
nl A cd 1 n n 

n n \ 

-( -1) n [q n < t -l/c> u ( t - 1 / G) - q n (t- -^ c ~ ) u(t-^ )J J ( 47 ) 

Hence, equation ( 41 ) reduces to 

A (t) = j 2 sin co L/ 2 c (AR cos [co (t-L/ 2 c) + y ] 

n A co n \n n n 

n n 1 

- (-l) n BQ cos [ co (t- lf ^ L ) -y ]) 

n n c 'n 7 

+ 2 sinix l(ar ’ cosine (t-L/ 2 c) + 5 7r] 

^ n V n n 

- (- 1 ) n BQ^ COS [ 0 n (t- — - -f^~) -i 7 T]) 

-AR" (e XnL -l)e" 9 n t 
n 

+ C 2 + R sin( co t + y ) + R' sin ( © + 5 n ) + R" e 9 n | 

L n n n n n n J 

- (-i) n D(l +Q n sin[co n (t-l/c) * 7 n 1 + Q^sin [©^t-i /c) - ]) J ( 48 ; 

where A, B, C and D have the following tabulated values 




A 

B 

C 

D 

0 < t < — 
c 

a 

0 

0 

1 

0 

L 1 


1 




— <t < — 
e e 

b 

0 

0 

0 

1 1 + L 

< t < 

c c 

c 

1 

0 

0 

1 

1 +L * 
< t 

c 

d 

1 

1 

0 

0 
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RESONANCE OF PLATE 


There are three conditions which may cause resonance of the plate to take 
place. 

1. Frequency resonance, = krc 

For the case of impulse load, the functions ip ^ and ip become very large 

if the damping factor is very small and for distributed load A^(t) , given by 
equation (41) , become infinite. The critical values of r, T=27rk/co n . versus 
the plate aspect ratio are plotted in Figure 6. 

2. Velocity resonance, /u =1 ore =co /A 

J n n n n 

When the frequency of the velocity waves produced by the traveling loads 
coincides with the natural frequency of the plate, the coefficients and be- 
come very large for a small damping factor, as shown in Appendix B for the 
case of impulse loads. For distributed loads, it can be seen from equations 
(38a) and (38b) that the inverse Laplace transforms will have a term 

L -1 F 2 2SC °2n 2 1 = tsinu; t . 

L (s z + u/) 2 J 

The critical values of c, c**= co /A > versus the plate aspect ratio are plotted 

n n 

also in Figure 6. The nondimens ional velocity of sound at standard atmosphere 

is c = c T/a s 0. 01 (a/h) for both steel and aluminum plates, 
s s 

3. Velocity -frequency resonance, ct = 2k 7: /A q 

Under this condition both g n (t) , given by equation (22b’), and A^(t) given 

by equation (41) , become unbounded. The critical values of r versus c are 
plotted in Figure 7. However, because of the smallness of the coefficient of g n (t), 

this condition should not impose serious response for the case of impulsive loads. 
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r = a/b 


FIGURE 6. c* AND t* VERSUS r 





EXAMPLES AND DISCUSSION OF RESULTS 


Examples on a plate subjected to a series of impulses: 

Example 1 . - Consider a steel plate with the following data 

a = 0. 6096m (24 in. ) , £ = 0. 1, v = 1/ 3 

E= 20. 67N/ m 2 ( 30 x 10 6 psi) , p = 7855 Kg/ m 3 ( 15. 24 slug/ ft 3 ) 

Consider also that velocity of the moving impulses is in the neighborhood of the 
speed of sound in a standard atmosphere. Now, examine how the stiffness of the 
plate affects the dynamic response by taking h = 0. 0254 mm (0. 01 in. ) and 
h = 0. 0803 mm (0. 03162 in. ) , respectively. 

Case I. h = 0. 0254 mm ( 0. 01 in. ) 

From the given data, calculate 

T = ~§- ^= 0.04186 T = 1 / f T= 0.09 c =cT/a' 23.41 
Ai D s s 

where c is the nondimensional velocity of sound. The values of ju = cX /co 
s n n n 

of this velocity for plate aspect ratio r = 1 to 4 are as follows 



i) 

r = 2 

r = 3 

r = 4 

n = 1 

68.7 

25.1 

11.84 

6.78 

n = 2 

55.7 

32.1 

17.56 

10.63 


The values of Aj and /3 are computed for the arbitrarily chosen velocities 
c = 15, 20 and 25. The results show that there is no appreciable difference 
for the three velocities because c is much greater than co /X^. Figure 8 shows 

and versus r. Note that for these velocities the approximate expression 
given by equation (34) is applicable. 

To calculate the upper bound of the maximum deflection and maximum 
bending stress, first read the value of A 1 and (3 from Figure 8 and and S 
from Figures 3 and 4, respectively; then use equations (26-1) and (30-1) . 

The results of w * and a* versus the plate aspect ratio are plotted in Figures 
9 and 10. 
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1 1. 5 2.0 2.5 3.0 3.5 4 . 0 

r 


0 


FIGURE 8. A t AND/3 VERSUS r OF CASE I, EXAMPLE 1 

Case II h = 0. 0803 mm ( 0. 03162 in. ) . Increase the thickness of the plate 
by \ITo times that of Case I, the values of T, c and are reduced to 1/ 10 

of the values calculated in Case I (t= 0.9). Plots of A t and /3 versus r for 

c = 1.5, 2.0 and 2.5 (c is 2.341 for this case) are shown in Figures 11 and 12 

s 

respectively. Figures 9 and 10 show the upper bound of the maximum deflection 
and maximum bending stress versus r. 

Example 2. - Considering again a clamped plate subjected to a series of 
impulses, assume that the nondimensional velocity of the impulses are small 

(t _ -) 

so that c « w /A and e 00 c } « 1 . By applying the approximation given 

n n 

by equation (35) , obtain A t and j3 for various values of r. The upper bounds 
of the maximum deflection and maximum bending stress are shown in Figure 13. 
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1 1.6 2.0 2.5 3.0 3.5 4.0 


FIGURE 13. w* AND a* VERSUS r FOR n«i 

Discussion of Results. - For a given frequency of the moving impulsive 
loads, the effects of the thickness of the plate, the plate aspect ratio, and the 
velocity of the loads on the dynamic responses are illustrated in Figures 8 to 
13. Some interesting results have been observed: 

1. If c» co /X (i.e. , with high velocity loads moving on a thin plate) , the 

n n 

dynamic response of the plate is almost independent of the velocity. 

2. If c« co n / A. , the dynamic response is directly proportional to 
the velocity. 

3. The three conditions of resonance with either n = 1 or 2 are: 
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(2) t = 2kn / co 

' ’ n 

(3) ^cX^t = k 7 r (this condition has no significant in reality) 

For a plate with large damping, condition (1) is a predominant factor; 
while condition (2) is more important for small damping. 

4. The peaks of the curves in Figure 11 correspond to the critical values 
of c and r given by Figure 6. 

5. The value of /3 (= A^Aj) is small for a plate with panel aspect ratio close 
to one. This is an indication that the two-term approximation used in solving 
equation (13) is adequate (more terms are required to cope with the higher modes 
if r is large) . 

6. Since the w* and <x* presented here are the upper bounds of the dynamic 
response; in reality, these values could be considerably higher than the actual. 
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APPENDIX A. FORMULATION OF a n (s) 


By definition of equation ( 8) , 



-sx/ c 
e 


[chAx - cosAx - a 


( shAx - sinAx ) ] dx 


-s/ c 
e 

shA - sinA 


j(- chA - AshA) (shA - sinA) + (fr shA + AchA) (chA-cosA)J 


+ (f - o' A) } / (s 2 /c 2 - A 2 ) 


- i e 


-s/c r s 


j(- ■§ cosA + AsinA) + a (“ sinA + AcosA)J + ^ - a xjy^sVcV A 2 ) 


Using the approximations 


n 


a =1, cosA =0, sinA = (-1) , and chA =shA , 

n n n n n 


obtain 


2c(cA ) 2 


a n^ (s + cA )(s 2 + c 2 A 2 ) 


n n 


, / s + cA 

. ..n -s/cf n 

- ( -1) ce l 0 *, 2 

\ s z + c'A 
x n 
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APPENDIX B. FORMULATION OFf(t) AND git) 


The terms contained in equation (19) are tabulated as follows: 


Laplace Transform Pairs and Coefficients of R( s) and Q( s) 


UL 

E4EH 

V z > 


Q k 

B 

l 

e- 0z 

1 


s + 0 

A 

e -1 

1+/u 2 -2ZM 

0 

I 

s + 

-£ on Z -Za>( z+t) 

e coscoh (z + t) - e coscorfz 

2m 2 R 0 BD 

-CE 


2(cosh£a>T - cosu^t) 


1 H 

-Zco z . -£c o (z+t) . 

e sinwH (z +t) - e sinatfz 



o 

■K WM 

2m 2 R 0 BF 

CG 


( s + £c o) 2 +co^ 

2(coshfu>r - cosco^r) 

o 

s 

-icsc|0r sin0 (z + ^t) 

-BD 

BE 

o 

s 2 + 0 2 

A 

0 

|csc 50 T cos© (z+ Jt) 

BE 

BD 

4 

s 2 + 0 2 


Definition of Constants 


b = [(i V ) 2 + 4 t,Vr 1 

c = 

[(1 V ) 2 + 8 £,V] ■ 

D = 1 - jU 2 + 2t,/u 

E = 

1 - M 2 “2^M 

II 

“E 

l 

1-^ 

1 

i M 

o: 

G = 

m(m 2 - i) +Um 2 +1 ) 


Notice that the subscripts "n" are omitted from all the symbols except s, z, 
and r, and that the constants B to G are approximations after neglecting the higher 
order terms of £ . 
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Further simplification by combining 4>i, *h and 4>z, respectively, 
results in 

f r (z) =R^i(z) +R 2 fc(z) (B-la) 

_ ( -Ccoz , . . „ , -£,co(z + t) . .,) 

= R|e cos[co ^(z+ r ) -£]-e cos (co ^z - £ ) | 

f (z) = Q^z) +Q 2 if>z( z) (B-lb) 

H 

f - Sco z - £co ( z+ t) . .1 

= Q|e cos [co ^ (z + r) - e] - e cos(co ^z - e) / 

g r (z) = R 3 ^ (z) + R 4 ^ 4 (z) =iN/2Bcscier sin [9(z+ %r)+ y ) (B-lc) 


g (z) = Q 3 ^3 (z) + Q 4 i/) 4 (z) = ^/2Bcsc^0 t cos|0(z+ £t) + y\ 

q L J 


(B-ld) 


where 

R 

Q 


jU 2 RnBN/p 2 + F 2 

cosh f CO T - COSCO .T 
d 


cj E 2 + G 2 


2 ( cosh f co t - cosco ,r ) 
d 


= tan -1 ( E/ D) . 


£ = tan -1 (F/D) 
e = tan -1 (G/ -E ) 


Substituting from equation (B-la) into equation (22c) and combining terms, 
we obtain 


A A 

f(t) 


A 

, a a _ t i- a n 

R 0 ^ 0 (t)+cje w cosL (t+T)-J 


-S^(t+r) . . 

- e cos (a>-t 

a 


where 

A 

c 


-<*>)) 


( B-2) 




( f > = tan' 


,2 , ^2 a -2£u(l-T-c)_ 9( _ 

Rsin| - (~l)“Qe 


( tT - r ) 


R z +Q 2 e — c , _ 2 (-i)“RQe - c ' C os ( £-e+ co it-%) 

1 


n^- £co( it - <r) . _2L. 

° sin( g+c - a) lt) 


„ . n - £co( lt - c ) / co . 

Rcos£ - (-1) Qe ^ cos( e + — - co it) 

c 
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Then, 

A A 


A A A A -too £ A 

f(t) = R 0 (t; + Fe cos(co d t -(f)- ip), 


where 

A 


AA 


F = CD, D = J” (e T -cosco r) 2 + sin 2 w,T , 


ip = tan 


-l 


sino) ,r 
d 


- Jc C T 

e - cosw.t 
d 


( B-3) 


( B-4) 


( B-5) 


It follows from equation (22d) that 


g(t) = | '^ 2 Bcsc| 0 t |sin [9 (t + ) +7 ] -(-l) n cos [^9 (^+ ^t+ tT - 3 )+y]j 

= 2^2Bcsc5 9 T^Gsin[ 0(^+ ir) +y + v]j . 


where 


A 

G 


= 4 


1 + (-l) n sin0(r - )| 2 +[cos0(t - ^ )J 2 

= J 2{ 1+ ( —1) n [ sin9rcos\ - cosGrsinA] } ^use cos 0, sinA n “ ( -1) 


= 2\/j (1 - cos 9 t) = 


2sin|9T 


rj = tan -1 


- (-l) n cos9(t- ~ ) 

V 

l-(-l) n sin0(0r- ~ — ) 

c 


= tan -1 


-sin9r 
1+ cos9t 


Hence, 

g$>= AT2Bsin [0(t + iT) + y+ rj] . 


(B-6) 


35 



APPENDIX C. FORMULATION OF EQUATION (40) 


In the following, reformulate the typical terms contained in equation (40) 
into simple form; 


1 . The unit step functions 


fi 


[u(t - kr) - u(t - Rt - L/c) ] =• 
£u (t - kT - i/c) - u( t - kr - -S]- 


0 

1 

0 


kr st <kr + L/c 
kr+L/c st s (k+ i) t 
kr+ i/c <t <k T +(L+ i)/c 
kr + (L+ i)/c s t < (k+ i) T +i/c 


2 . The exponential functions 
K 

Y exp |"-c\(t - kr)l u(t - kr) - exp 
k=0 L J 



kr - L/c)J u(t - kr- L/c) 


-c\t(K +1) 

, XL 1-e -c\t n 

/ - (e -i) : e 

1 - e 

l - (e AL -l) - 1 " e 


-cXr 

-cAtK 


1 - e' cXT 


L/cgt 0 S r 
-cMt,+r) e +e -cxt, 0£ , < L/0 


where t 0 = t - Krwith Kr<t< (K+ 1 )t . 

3. The sine functions 


Making use of the identity, 

K 

Y sink© =sin|(K + i)© sin|K0 cscje, 
k=0 

we obtain 

Y sin[ 00 (t-kT) +«] u(t - kr) - sin[co (t-kr-L/c) +a] u(t-kr- L/c) 
k=0 
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^ 2sin-~~ | sin ^Ko) t csc^w t cos |a) (t 0 + ^pr) + 4>J + cos[w(t 0 + Kr) + <£]J 

if Kt + L/c < t < (K+ l)r 

^ gsin^— 1 jsin^ (K - 1 )o>tcsc gcoTCOs [co (t 0 + |Kt )+<£]+ cos[co (t 0 + Kt) + <pjJ 

+ sin (oj $ 0 + a ) , 

if Kt < t < Kr + L/c, 
where <{>=a - uL/2c. 
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